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Abstract 

We study cosmological perturbations in the context of an interacting dark energy model, in 
which the cosmological term decays linearly with the Hubble parameter, with concomitant matter 
production. A previous joint analysis of the redshift-distance relation for type la supernovas, 
barionic acoustic oscillations, and the position of the first peak in the anisotropy spectrum of the 
cosmic microwave background has led to acceptable values for the cosmological parameters. Here 
we present our analysis of small perturbations, under the assumption that the cosmological term, 
and therefore the matter production, are strictly homogeneous. Such a homogeneous production 
tends to dilute the matter contrast, leading to a late-time suppression in the power spectrum. 
Nevertheless, an excellent agreement with the observational data can be achieved by using a higher 
matter density as compared to the concordance value previously obtained. This may indicate that 
our hypothesis of homogeneous matter production must be relaxed by allowing perturbations in 
the interacting cosmological term. 
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I. INTRODUCTION 



The cosmological constant problem has acquired a renewed importance since several inde- 
pendent observations have been pointing to the presence of a negative pressure component in 



the cosmic fluid 



From the point of view of quantum fleld theories, the natural candidate 
for such a dark energy is the quantum vacuum. Since, at the macroscopic level, it has the 
symmetry of the background, its energy-momentum tensor has the form = Agf^, where A 
is a scalar function of coordinates. This leads, in the case of an isotropic and homogeneous 
space-time and co-moving observers, to the equation of state ptv = — Pa = —A, where A may 
be, in general, a function of time. In the case of a constant A, the vacuum contribution 
plays the role of a cosmological constant in Einstein's equations. 

However, the estimation of the vacuum energy density by quantum fleld theories in the flat 
space-time leads, after some regularization procedure, to a very huge result when compared 
to the observed value. A possible way out of this difficult is to argue that such a result is 
valid only in a flat background, in which the very Einstein equations predict a null total 
energy-momentum tensor. Therefore, the huge vacuum density should be canceled by a 
bare cosmological constant, like in a renormalization process. Now, if we could obtain the 
vacuum density in the FLRW space-time, after the subtraction of the Minkowskian result it 
would remain an effective time- dependent A term, which decreases with the expansion. 

The idea of a time-dependent cosmological term has found different phenomenological 
implementations 2[], being a subject of renewed interest in recent years (3|, y, Isj. A general 
feature of all those approaches is the production of matter, concomitant with the vacuum 
decay in order to assure the covariant conservation of the total energy |^. Indeed, in the 
FLRW space-time, the Bianchi identities lead to the conservation equation 

PT + 3i/(pr + Pr) = 0, (1) 

where pT and pt stand for the total energy density and pressure, respectively, and H = a/a 
is the Hubble parameter. By writing px = Pm + ^ and = Pm — (where pm and Pm are 
the energy density and pressure of matter), the above equation reduces to 

pm + SH{pm + Pm) = "A, (2) 
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which shows that, in the case of a varying A, matter is not independently conserved^. 

An important point to be clarified in this kind of model is the homogeneity of matter 
production. Of course, in a strictly homogeneous space-time the production is homogeneous, 
since pm and A depends only on time. But, in the presence of density perturbations, is the 
new matter produced homogeneously, or just where matter already exists Q|? In the case 
of a homogeneous production, the new matter tends to dilute the density perturbations. 



leading to a suppression of the density contrast. In some mode 
enough to impose very restrictive observational limits to them 



^ this suppression is strong 
9|. 



In this paper we will analyze the evolution of density perturbations in a particular. 



spatially flat, cosmological model with vacuum decay |10l . It can be based on a phe- 
nomenological prescription for the variation of A with time jl2| . given by A ^ (if +m)'^ — m^, 
where m is a characteristic energy that can be identifled with the scale of the QCD vacuum 



condensation, the latest cosmo 



ogical vacuum transition. Although it can be corroborated 



by holographic arguments [12, ll3| , based on the thermodynamics of de Sitter space-times, 
here we will take it just as a phenomenological ansatz. In the limit of very early times, we 
have A ^ H"^, which provides a non-singular inflationary solution [l^ . 

In the opposite limit of large times we have A = cxH, with a ~ m^. This scaling law 
for the vacuum density was also suggested in J^], on the basis of different arguments. It 
leads to a cosmological scenario in qualitative agreement with the standard one jl^, with 
an initial radiation era followed by a long phase dominated by dust. This dust phase 
tends asymptotically to a de Sitter universe, with the deceleration/acceleration transition 
occurring some time before the present epoch. On the other hand, a quantitative analysis 
has shown a good accordance with supernova observations, leading to age and matter density 
parameters inside the limits imposed by other independent observations 

Since the radiation phase we obtain is indistinguishable from the standard one, our anal- 
ysis will be initially focused on the evolution of density perturbations of non-relativistic 



^ Properly speaking, we should also consider the pressure and energy associated to the very process of 
matter production, that is, the energy-momentum tensor of the interaction between matter and vacuum. 
In this sense, decaying vacuum models do not differ essentially from interacting dark energy models 0, 
with the scalar function A replaced by a scalar field interacting with matter. Nevertheless, if the vacuum 
decays into non-relativistic particles, as we will consider here, the interaction term can be neglected, and 
the above decomposition may be considered a good approximation. 
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matter in the dust-dominated phase, considering wavelengths inside the horizon^. In this 
way, it will be possible to make use of a generalization of the Newtonian linear treatment of 
the problem, which includes the effects of matter production We will show that, even in 
the case of a homogeneous vacuum decay, the contrast suppression is important only for late 
times, not affecting the process of galaxy formation. On the other hand, it dominates for 
future times, and we will discuss how this behavior can possibly alleviate another problem 
related to the cosmological term: the cosmic coincidence problem. 

Subsequently, a relativistic analysis will be performed, in order to construct the matter 
power spectrum. Again, the hypothesis of homogeneous matter production will be used, 
leading as well to a consequent power suppression. A second interesting difference as com- 
pared to the ACDM model is a shift of the spectrum turnover to the left, that is, to smaller 
wavenumbers. The late-time suppression is not very sensitive to the value used for the mat- 
ter density, a feature that can already be noted from the Newtonian analysis. On the other 
hand, the correction of the turnover position, by taking a higher matter density, displaces 
all the spectrum to the right, compensating the late-time power suppression. In this way 
we can obtain an excellent fit of data, but with a higher matter density in comparison with 
the standard case. 

The article is organized as follows. In next section we review the main features of our 
interacting model. In Section III we perform the Newtonian analysis of evolution of density 
perturbations in the matter era. In Section IV the matter power spectrum is constructed, 
on the basis of a simplified relativistic calculation. In Section V the reader can find our 
concluding remarks. 

II. THE MODEL 

The Friedmann equations in the spatially fiat case are given by ([T]) and pt = 3H^. Let 
us take px = Pm + Pt = Pm — and Pm = {l — ^)Pm, with constant 7. Let us also take 

^ As already commented, we will assume that the vacuum is decaying into non-relativistic particles, in 
order to avoid any conflict with CMB observations and with the observed coldness of dark matter. We 
will also suppose that only dark matter is produced, since the baryon content is well constrained by 
nucleosynthesis. Evidently, these assumptions cannot be verified without a microscopic theory of the 
vacuum-matter interaction. 
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the ansatz A = aH^ with a constant and positive. We obtain the evolution equation 

2H + 37iJ2 - a^H = 0. (3) 



The solution, for p^, -f^ > 0, is given by |lO[ | 



a = C[exp((T7t/2) - 1]^ , (4) 

where a is the scale factor, C is an integration constant, and a second one was taken equal 
to zero in order to have a = for t = 0. 

In the radiation phase, taking 7 = 4/3 and the limit of early times {at << 1), we have 



a ^ ^2CVt/3. (5) 

This is the same scaling law we obtain in the standard case, leading to H ^ l/2t. In the 
same limit we then have pm = Pt — ^ = 3if^ — aH ^ 3if^ = Pt- By using ^ we then 
obtain 

PT-P.-^-4^, (6) 

i.e., the same variation law for radiation one obtains in the standard model, which shows 
that, during the radiation era, both the cosmological term and the matter production can 
be dismissed. 

On the other hand, in the matter era we obtain, by doing 7 = 1, 

a = C[exp(aV2) - 1]^ . (7) 

Taking again the limit of early times, we have 

a ^ C{atl2f'\ (8) 

as in the Einstein-de Sitter solution. It is also easy to see that, in the opposite limit t — > oo, 
([7]) tends to the de Sitter solution. 

With the help of ([7]), and by using A = aH and pm = — aH, it is straightforward to 
derive the matter and vacuum densities as functions of the scale factor. One has 

^"^ ~ 3a3 + 3a3/2 ' ^9) 



In these expressions, the first terms give the standard scahng of matter (baryons included) 
and vacuum densities, being dominant in the limits of early and very late times, respectively. 
The second ones are owing to the process of matter production, being important at an 
intermediate time scale. 

With and fllOp we obtain, for the total energy density and pressure^, 

a J 



a 

Pt = 



- +1 



(11) 



Pt = -^YpV"- (12) 

From (I7j) we can also derive the Hubble parameter as a function of time in the matter 
era. It is given by 

H = ^7^-^- (13) 

1 - exp(-(Tt/2) ^ ^ 

With this expression, and by using ([7]) and ([9]), it is not difficult to obtain the present age 
of the universe, given, in terms of the age parameter, by 

Hoto = -^^^^ (14) 

3("mO ~ 1) 

where Qmo = Pmo/S-^o relative matter density at present. 

Finally, with the help of (I7j) and ( |T3l) . we can express as a function of the redshift 
z = ao/a — 1, which leads to 



Hiz) = m 



l-Q^^ + Q^,(z + lf/'] . (15) 



With this function we have analysed the redshift-distance relation for type la supernovas 



ll| , obtaining data fits as good as with the fiat ACDM model. With the Supernova Legacy 
Survey (SNLS) [18] - the most confident survey we have so far - the best fit is given by 
h = 0.70 ± 0.02 and n^o = 0.32 ± 0.05 (with 2a), with a reduced x-square Xr = 1-01 
(here, h = i/o/(100km/s.Mpc)). On the other hand, a joint analysis of the Legacy Survey, 
baryonic acoustic oscillations and the position of the first peak of CMB anisotropies has led 
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These are the same expressions we obtain for a generahzed Chaplygin gas (characterized by the equation of 



state Pch ~ '^^/Pch if choose a — —1/2 and A = yV^/S (see for a detailed discussion about 

this and other curious equivalences between dark energy models). Note, however, that the oscillations in 
the evolution of density perturbations characteristic of a Chaplygin gas are not present in our case, 
as we will see below. 
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;o the concordance values h = 0.69 ± 0.01 and Qmo = 0.36 ± 0.01 (with 2a), with Xr = 1-01 
lol . With these results one can obtain, from f|T^ . a universe age to ~ 15.0 Gyr, inside the 
interval allowed by age estimations of globular clusters [20 1. 



III. NEWTONIAN EVOLUTION OF DENSITY PERTURBATIONS 



The Newtonian equation for the evolution of density perturbations in a pressureless fluid 



can be generalized in order to account for matter production 
it is given by 



2 pm Ot \pm 



In this generalized form, 



(5 = 0. (16) 



Here, 5 = 6pm/ pm is the density contrast of the pressureless matter, and \1/ is the source of 
matter production, defined as 

p^ + 3H{prn+Pm) = '^- (17) 

In the case of a constant A, \1/ = 0, and f|T6|) reduces to the usual non-relativistic equation 
for the linear evolution of the contrast. In our case, on the other hand, = —A = —aH, as 
can be seen from ([2]). 

Equation (fT6l) is derived on the basis of two main assumptions [14]. The first one is 
that the produced particles have negligible velocities as measured by observers co-moving 
with the cosmic fluid. This is a reasonable hypothesis, since we are dealing with a non- 
relativistic phase of universe expansion, when H (and so A) varies slowly enough. The 
second assumption is that the vacuum component A is strictly homogeneous, which means 
that matter production is homogeneous as well. This stronger hypothesis is totally ad hoc 
at the present stage of the model development and, as we will see, leads to a suppression of 
the contrast at large times. 

In order to solve f fT6l) for our case, it is convenient to introduce the new variable 



X = exp{-at/2). (18) 

After calculating pm, H and ^ as functions of x with the help of (JTj), ([9]) and (|T3i) . equation 
(fT6|) takes the form 

?>x^{x - If 5" + Ax{x - 1)5' - 2(3a; - 2)5 = 0, (19) 



where the prime means derivative with respect to x. It is possible to show that, in the hmit 
of early times, it reduces to the evolution equation for the contrast in the Einstein-de Sitter 
model, as should be. 

The general solution of f[T^ can be written as 



X — 1 



1/3, 2/3) + 1)2/3 

3 



(20) 



where Ci and C2 are integration constants to be determined by initial conditions, and 
j3{x, a, h) is the incomplete beta function, defined as 



/3(x,a,6)= / y'^-' {l~yf-' dy. (21) 







This (3 function can be expanded in a Laurent series around a; = 1, leading to 

/3(x, 1/3, 2/3) ^ /3(1, 1/3, 2/3) - ^(x - 1)^/^ (22) 

In this way, with the help of fl22|) and (fT8|) we can expand (120!) around t = 0, obtaining 

d^^ + D^e/^ (23) 

which is precisely the general solution obtained in the Einstein-de Sitter model, as expected. 
The new arbitrary constants are given by 



a 



Ci + ^/3(l, 1/3, 2/3)^2 



(24) 



D2 = ^-J^C2. (25) 

If, in the early time approximation (!23l) . we want to retain just the growing mode, pro- 
portional to t^/s^ must choose Di = 0. Then, our general solution fl20l) reduces to 

± ^ 2x[/j(l,l/3,2/3)-/?(x,l/3,2/3)] _ x^/^ 

C2 3(1 -x) (1 -x)V3- ^ '^^ 

The above solution can be expressed as functions of t or a, with the help of (1181) and ([7]). It 
can also be expressed as a function of the redshift, by using the relation 



which can be derived with the help of ([71), ([9]) and (U31) . 

Figures 1 and 2 show the density contrast (l26l) as a function of a and z, respectively. We 
have taken ao = 1, and used for the matter density parameter the best-fit value we have 
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FIG. 1: The density contrast as a function of the scale factor 



obtained from the SNLS analysis Uj, flmo = 0.32. The integration constant C2 was chosen 
so that for the time of last scattering {z ^ 1100) one has 6 ~ 10^^, as imposed by anisotropy 
observations of the cosmic microwave background j2l|. For the sake of comparison, we have 
also plotted the evolution of the density contrast in the Einstein-de Sitter solution and in 
the spatially flat ACDM model with fi„o = 0.27. 

In our case the density contrast grows monotonically with time until z ^ 0.6, after which 
it decreases monotonically, tending to zero in the limit t — 00. The consequences of such 
a suppression at large times will be discussed in our Conclusions, where a possible relation 
with the cosmic coincidence problem will be outlined. The important point here is that the 



. 05 




FIG. 2: The density contrast as a function of the redshift 
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evolution of S in our case is indistinguishable from its behavior in the ACDM case until 
z ~ 5, that is, along the entire era of galaxy formation. On the other hand, the late-time 
suppression leads to a present contrast approximately 1/3 of the standard one, a difference 
that will be manifest in the power spectrum, as we will see now. 



IV. THE POWER SPECTRUM 



The shape of the spectrum depends on several parameters. But one of the most important 
is given by the moment of equilibrium between radiation and matter, Qji = Q^, where Qji 
and Qrn are the respective density parameters (relative to the critical density). In the ACDM 
model, we have 

nn = ^ = nno{l + z)\ (28) 



n^ = ^-^ = nmo{i + zf, (29) 



a3 



where Qjiq and Qmo are the density parameters for radiation and matter today. The redshift 
at equilibrium is then given by 



Following 



1 + .e. = (30) 



22], we fix, for the ACDM model, 

fi^o/i^ = 0.127, riRoh^ = A.l X 10-\ (31) 



This implies 

1 + Zeg = 3097. (32) 

Remark that this value, as a matter of fact, is independent of h. 

Now, we can analyse the moment the perturbations enter in the horizon. This is obtained 
by inspecting the perturbed equations. In general, it can be written as 

In this equation, 6 is the density contrast, and Vg = ^ represents the sound velocity in 
unities of c (the velocity of light). The presence of a first derivative term is related to the 
friction due to the expansion of the universe, while the two last terms describe the interplay 
between the pressure, that avoids the collapse, and the gravitational attraction, that drives 
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the collapse. When the first of these terms dominates, the perturbation does not grow; when 
the second one dominates, the perturbation increases. Ignoring numerical factors of order 
of unities, related to the sound velocity, equation of state etc, the condition that separates 
both regimes is 



, a c ca 

k=-r^ dH = -fj = —- 
clh Ha 

In this expression, dn is the Hubble radius. Of course, this is just an estimation. 
For the ACDM model, we have 

-1/2 



(34) 



dn = — n^oil + zY + Qnoil + zf + ^ao 

-no 

where I^ao is the density parameter for the cosmological term today. Hence, we have 



[(1 + z)kiH^f = a„o(i + zf + fi^o(i + zf + n 



AO; 



(35) 



(36) 



where Iho is the Hubble's radius today, Iho = 3000 /i"^Mpc. In general, for large values of 
z the term ^ao can be ignored. In doing so, and using the expression above for z = Zeq-, we 
find the formula (7.39) of reference 



23| 



2 n 



mO 



Using, besides the values of VL^q and VLrq already quoted, also /i = 0.7, we obtain 



fceg = 0.013. 



(37) 



(38) 



We notice that, using the BBKS transfer function for the ACDM model 2J], the turning 
point is also located at = 0.013. 

Now, the observations cover scales from kminh~^ = 0.010 until kmaxh~^ = 0.185. Using 



the parameters above, we find that these modes entered in the Hubble horizon at 



kmin ^ Zi 2077, kfnax 



Z2 = 59143. 



(39) 



That is, essentially, all modes entered in the radiation dominate era. 

Turning to the present interacting model, the main modifications are the following: 

1. The expression governing the moment the modes enter in the Hubble horizon is given 

by 



klm (1 + z) 



1 



+ nnoil + z) 



(40) 
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with flrno + ^Ao ~ 1- This is an approximate expression obtained from (fTSi) by adding 
a conserved radiation density to the Friedmann equation 3H^ = px-^ 



2. An inspection of (fT5l) for high z, when A and the matter production are dismissable, 
shows that ^lm.{z) = Q'^(^{l+zY . In other words, we have the same scahng of conserved 
matter as in the standard model, but with an extra factor Qmo- This is owing to the 
matter production between t{z) and to: in order to have the same matter density 
today, we need a smaller density at high redshifts. As a consequence, the redshift of 
equilibrium between matter and radiation is now given by Zeg = Q'^q/Qrq, while for 
the correspondent wave number we obtain, instead of fl37|) . 



V "fio I'm 

Note the extra factor Qmo as compared to the corresponding ACDM expression. As 
this factor is smaller than unity, this means that the turnover of the spectrum is moved 
to the left, that is, to smaller fc's as compared to the standard model. 

3. The matter density parameter and the Hubble parameter are not the same as before. 
In the subsequent analysis we will use Qmo = 0.32 and h = 0.7 (the type la supernovas 



best fitting 
Now, the results are the following: 

1. The equilibrium occurs at = 2263, which implies k^g = 0.007; 

2. The mode kmin enters in the Hubble horizon at zi = 3469, while the mode kmax at 
Z2 = 81404. 

As already noticed, the results indicate that the spectrum is displaced to the left, implying 
that there is a power suppression with respect to the ACDM model. Moreover, there is, as 
we have seen in the previous section, an additional power suppression during the matter 



^ Note that the inclusion of conserved radiation changes the dynamics, and, consequently, the production 
of matter, A(z) and Pmiz) also change. Therefore, the exact generalization of (jlSp requires a reanalysis 
of the dynamics. Nevertheless, as ^ro ~ 10"** << 1, when the vacuum and the matter production begin 
to have importance, the radiation is negligible, and vice-versa. In this way, (|40|) can be considered a very 
good approximation. Indeed, a numerical analysis in the range < z < 10^ has shown that the difference 
between (|40|) and the exact H{z) is as small as 0.01%. 
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dominated phase. Hence, essentially, we must expect that the power spectrum displays, in 
what concerns matter agglomeration, an expressive power suppression in comparison with 
the ACDM model. 

However, we can displace the spectrum to the right, instead of displace it to the left, if 
the values of flmo and/or h are increased. For example, for = 0.48 and h = 0.73, the keq 
occurs at 0.016, with Zeq = 5094. Moreover, kmin enters in the Hubble horizon at zi = 2589 
and kmax at Z2 = 80020. The substantial displacement to the right of keq compensates 
the smaller growing of perturbations during the matter dominated phase. So, the general 
features of the power spectrum are reproduced for larger values of Qmo as compared to the 
ACDM model. 

A precise derivation of the spectrum is a very though calculation, since the Einstein- 
Boltzmann coupled system must be considered. A complete analysis for the ACDM model 
leads to the so-called BBKS transfer function [2^, which gives the spectrum today as func- 
tion of a given primordial spectrum. For the scale invariant spectrum, favored by the 
primordial inflationary scenario, the BBKS transfer function is given by 

PUk) = \SUk)\' = AT(A:)4^fc, (42) 

where ^4 is a normalization of the spectrum (which can be fixed by the spectrum of anisotropy 
of the cosmic microwave background radiation), T{k) is given by 

ln(l + 2.34g) 



T{k) 



2Mq 



1 + 3.89g + (16.1g)^ + (5.64g)-' + (6.71g) 

k ^ n 



(43) 



q = J^^PC~\ r = n,,^ohe ^^mo , (44) 



and where ^mo, ^dmo, and Qt are, respectively, the present density parameters of 
pressureless (baryonic + dark) matter, dark matter, baryons and the total energy. The 
function g{Q) is defined by 
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n7 - ^]Ao+( 1 



AO 



2 / V 70 



-1 



(45) 



The transfer function defined above represents the fitting of the complete numerical evalu- 
ation. 

A simplified version of the transfer function, which keeps all its essential features, can be 
obtained by integrating the perturbed equations for the coupled system containing radiation 
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and pressureless matter, from a very high redshift until today |25l. l26l|. The starting point is 
given by the Einstein equations and the conservation law for the energy-momentum tensor: 

R,, = ^-nGY.[TU-\9,ur], (46) 

Tr,, = 0, (47) 

where the indice denotes the ith fluid component. One of them will be radiation. The other 
one will be the pressureless matter in the ACDM case, or the vacuum-matter interacting 
fluids in our case (remember that, in our case, the pressureless matter is not independently 
conserved, since it interacts with vacuum). Introducing the perturbations, g^y = (7°^ + h^^^, 
Pi = p'i + Sp, Pi = p'i + Spi, with ((7°^, p°, p°) being the background solutions, and imposing 
the synchronous coordinate condition /i^o = 0, we end up with the following set of coupled 
equations: 

h + 2-h = p^6m + 2pRSn, (48) 
a 

Sm - —Sm = (49) 
Pm ^ 

(50, 

where h = hkk/o?-, and 5r are the density contrast for matter and radiation respectively, 
V is connected with the peculiar velocities of the perturbed radiative fluid, and in the ACDM 
case A is, evidently, zero. 

We now eliminate the variable h using fl49l) . divide all the expressions by Hq and rewrite 
the resulting equations in terms of the redshift which becomes the new dynamical variable. 
In the ACDM case the system of equations is reduced to 

^'R-\{^ + ^'m]=^. (53) 



where the primes indicate derivative with respect to the redshift z. The background functions 

f[z\ and g[z\ are given by 



a 



2 



a? 



f[z] = ^ = a^o(l + zf + ^Roil + zY + r^Ao, (55) 
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9[z] = ^ = -]pmo{l + zf - + zf + Qao. (56) 

Integrating, for example, from z = 10^ (when the initial spectrum is supposed to be scale 
invariant, i.e., 6m, Sr oc \/k) until today, z = 0, we can reproduce the BBKS transfer function 
with about 10% of precision. 

We can perform the same calculation for the present model, finding the following set of 
perturbed equations: 

,n r^A , 9i[z] 1 \ , (gi[z]n'^ 1 ni n'^n'^^ 
fi[z]{i + z)^i "'^\Mz]nm{i + zy ni / 

3 1 



+ 2nR6Ry (57) 



v' = -i'-^) (59) 

il + z)JMz\ 



In these equations, we use the following definitions: 

/iW = ^'= fo ^o\n^ {nAo + nmoil + zf'V 

(^AO + i^mo)(l + Zy [ } 

+ (l + zfnno, (60) 
9M = d = -^^±^f[[z], (61) 

^AO^mO f^mO/^ \3 i /'i i \- 



^m{z) = 1^(1 + ZY + (1 + .)n, (62) 

"AO + i 'mO "AO J 

^^aW = {l + ^(l + .)H. (63) 

i 'AO + i 'mO l- "AO -I 

In figures 3 and 4 we display the results for the exact transfer function for the ACDM 
model (blue), the corresponding numerical approximation (red) and the approximative trans- 
fer function for the present model (violet). The observational data come from the 2dFGRS 



galaxy survey program 



27|. In the case of the interaction model we used, in Figure 3, 



^mo = 0.36, the concordance value obtained from the joint analysis of type la supernovas, 
BAO and CMB [l^l. In Figure 4, on the other hand, we have used VtmQ = 0.48. We see that 
in the first case there is a substantial suppression of power, while in the second case, where 
the dark matter parameter has been increased, the agreement is excellent. 

Hence, concerning the matter power spectra, the interacting model with homogeneous 
matter production requires an almost double quantity of dark matter with respect to the 
ACDM model. 
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FIG. 3: The matter power spectra as given by the BBKS transfer function (blue), the approximative 
numerical analysis used here for ACDM (red) and for the interacting model (violet). The data come from 
the 2dFGRS galaxy survey program It has been used flmo = 0.36 for the interacting model. 

V. CONCLUSIONS 

In spite of the physical plausibility of a time dependent cosmological term, a complete 
theoretical development of this idea, including the microscopic details of the vacuum-matter 
interaction, is still lacking. On the other hand, macroscopic approaches depend on some 
phenomenological hypothesis, leading some times to diverse prescriptions for the vacuum 
decay^. For this reason, a careful comparison with current observations is very important, 
playing the role of corroborating or ruling out the different rnodels. 

We have already analysed the supernova observations [11], obtaining good fits and cos- 
mological parameters in accordance with other independent tests, as the age of globular 



clusters and dynamical limits to the matter density |28|]. Other precise tests, as the position 



^ For example, the linear dependence of A on the Hubble parameter we use here contrasts with the quadratic 
dependence used in reference In that work, the quadratic dependence is due to the computation of 
quantum effects of matter field in a cosmological background, which leads to a running cosmological term. 
The authors also used the matter power spectrum data to constrain the fundamental parameters of the 
quantum model. 
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FIG. 4: The matter power spectra as given by the BBKS transfer function (blue), the approximative 
numerical analysis used here for ACDM (red) and for the interacting model (violet). The data come from 
the 2dFGRS galaxy survey program {27I. It has been used flmo = 0.48 for the interacting model. 

of the first acoustic peak of the cosmic microwave backErround and the baryonic acoustic 
oscillations have also been performed [191], showing a good concordance when jointed to the 
supernova analysis. 

In the present paper we have studied the evolution of matter density perturbations, in 
particular the contrast suppression associated to the process of matter production. We have 
shown that, even in the case of a homogeneous production, the evolution of the contrast is 
the same as in the standard recipe along the entire era of galaxy formation, diverging from 
the later only for z < 5. 

On the other hand, the suppression would be dominant for future times, and this may have 
an interesting relation with another problem related to the cosmological term, namely the 
approximate coincidence between the present densities of matter and dark energy. Indeed, 
we can see from figures 1 and 2 that the matter contrast has its maximum just before 
today, when matter and vacuum give similar contributions to the total density. The largest 
structures formed until now tend to disaggregate in the future, and their existence then 
coincides with the time of approximate equality between the matter and vacuum densities. 



17 



This could alleviate the cosmic coincidence problem, if galaxies also follow such a process. 
However, we should remember that galaxies have left the linear regime of growth a long time 
ago, and that now their evolution is non-linear, driven essentially by their self-gravitation. 
Therefore, an explanation of the cosmic coincidence in the terms above will depend on a 
non-linear study of density perturbations in the context of the present model. Only such 
an investigation would tell us whether the contrast suppression described here can affect 
smaller structures like galaxies. 

Its also important to have in mind that the homogeneity of the matter production, implicit 
in the derivation of solution ( l26l) and in our simplified relativistic treatment, is just an ad 
hoc hypothesis, to be verified from both the theoretical and observational viewpoints. For a 
constant, non-interacting vacuum term it is certainly true, but not necessarily in the present 
case. Any inhomogeneity of the vacuum density around matter distributions may lead to 
an inhomogeneous production, reducing in this way the contrast suppression. This would 
allow us to fit the observed power spectrum with a smaller matter density, closer to the 
concordance value obtained in |l9| . Whether the matter contrast will still have a maximum 
around the present time, with the discussed implications for the coincidence problem, is a 
matter of investigation. A relativistic study of this case, that is, with 5A 7^ 0, is already in 
progress. 
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